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Abstract 

We consider a class of supersymmetric models containing baryon number violating processes 
such as observable n-n oscillations that are mediated by color triplet diquark fields. For plausible 
values of the diquark-quark couplings, the scalar diquark with mass between a few hundred GeV 
and one TeV or so can be produced in the s-channel at the LHC and detected through its decay 
into a top quark and a hadronic jet. 
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1 Introduction 



Being a proton-proton (pp) collider the LHC is an ideal hunting ground for new elementary 
particles which carry non-zero baryon number. Such particles often arise as vectorlike pairs 
in a large class of supersymmetric (SUSY) models, including some which are obtained from 
compactification of superstring theories [H El E]- 

The TeV scale vectorlike particles have various indirect phenomenological applications. 
They can help to substantially increase the mass of the lightest CP-even Higgs boson of mini- 
mal SUSY standard model (MSSM) [U E|. The new Higgs coupling to these particles can also 
modify the first order electroweak phase transition, and thereby make the scalar top quark mass 
bound milder in the electroweak baryogenesis scenario [6] . The leptonic vectorlike particles can 
generate sizable contribution to the tau and tau neutrino electric dipole moment which can be 
tested in future experiment [7]. 

The low scale vectorlike particle can also induce observable neutron-antineutron (n-n) os- 
cillations [3 [9], in which the so-called diquark couplings play a crucial role. The vectorlike field 
which couples to two anti-quarks have baryon number 2/3 (namely twice the baryon number of 
quark), and thus we refer to such a colored vectorlike field as "diquark". The coupling between 
diquarks and a pair of quarks can make the creation and decay modes of the diquarks very 
interesting in view of collider physics (TU]. If the diquark fields are pair produced at the LHC 
through QCD processes [HI [12], they decay into quarks via the diquark interaction. This in- 
teraction also plays a key role in the resonant production of the diquark field from pp collision, 
and if the cross section is large enough to distinguish the peak from the standard model (SM) 
background, one can detect the diquark field and determine its mass [13]. 

There are two possibilities for the representation of the diquark fields: (1) color triplet, and 
(2) color sextet. The color sextet diquark fields can be embedded in the Higgs multiplet required 
to break the gauge symmetry 577(4) c x SU{2) L x SU(2) R [H] down to SM, and the mass of 
the diquark can remain around a TeV scale if there exists an accidental global symmetry in the 
Higgs potential [8]. Such sextet diquarks do not have a renormalizable coupling with leptons, 
and their existence at low energy is not problematic. The color triplet diquarks are also well 
motivated since they can be accommodated into matter representations in unified models [2]. 
However, they can couple to both quarks and leptons, and the most general couplings can cause 
serious problems with rapid nucleon decay if the diquark mass is around the TeV scale. This 
difficulty is remedied by employing a U(l) symmetry to prohibit the couplings of the color 
triplets with leptons [9]. In this case, there are no dangerous nucleon decay operators, but the 
£7(1) symmetry can allow baryon number violating processes [9] like n-n oscillations. 

In this paper, motivated by the physics of n-n oscillations, we study the resonant production 
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and decay of the triplet diquarks at the LHC. The attractive feature of gauge coupling unifica- 
tion in MSSM is easily maintained |3j in the case of triplet diquarks accompanied by electroweak 
doublets. Clearly, the diquark couplings with the up and down quarks should be large enough 
to detect its resonant production. However, some of the couplings (especially for the first and 
second generations) are constrained by the experimental data on meson-antimeson (K-K, B-B 
and D°-D°) mixings, and thus we will first discuss the experimentally allowed parameter range 
of the diquark couplings. In the case of sextet diquark, D°-D° mixing is generated even at tree- 
level, while for the triplet diquarks, the mixings are generated by box or penguin diagrams, 
quite similar to the conventional analysis in MSSM. We will study whether one can see the res- 
onant peak from the SM background in the phenomenologically allowed region of the diquark 
coupling. If the triplet diquarks couple with quarks for both first and third generations, one 
can observe resonant decay events for both hadronic dijet and single top quark plus one jet. 
Note that the triplet diquark can couple to both left- and right-handed quarks, while the sextet 
diquarks can couple only with the right-handed quarks. It is therefore important to observe 
top quark helicity to distinguish these two cases of diquark production from each other. 

This paper is organized as follows: In section 2, we introduce following j9] a U(l) symmetry 
which forbids rapid nucleon decay processes but allows n-n oscillations. In section 3, the 
phenomenological constraints of the diquark couplings are considered. The diquark can couple 
with the quarks of all generations. In section 4, we study the collider phenomenology of the 
triplet diquarks. We first study the mass bounds on diquarks arising from the Tevatron results, 
and then discuss the resonant production and decays of the diquarks at the LHC. Section 5 
summarizes our conclusions. 

2 A model with triplet diquark couplings 

We introduce the following vectorlike fields (more precisely chiral superfield) with masses in 
the few hundred GeV to TeV range: 

D : (3, l)_i + D : (3, l)i + L : (1, 2)_i + L : (1, 2) i . (1) 

Because D + L has the same particle content as the 5-dimensional 577(5) multiplet, gauge 
coupling unification in MSSM is preserved. The vectorlike fields can couple with the quark and 
lepton fields as follows: 

Dqq + Du c e c + Dqt + Du c d c + L£e c + Lqd c + Lqu c . (2) 

Here we use the standard notation (q, u c , d c , £, e c ) for the quarks and leptons. In general, the 
couplings in Eq. (j2J) contradict the current experimental data [15] unless the coupling coefficients 
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Table 1: Particle charge assignments under (nB + mL)/2 symmetry. The numbers n, m, Xd 
and X L are all integers. 



are less than O(10 -13 ). Indeed, the exchange of the scalar component of D (D) chiral multiplet 
can generate the dimension-six proton decay operator qqu c *e c ^ (q£u c 'd c *), which leads to an 
unacceptable proton decay rate. 

The low scale vectorlike multiplets (5 + 5) can, in general, generate neutrino masses through 
loop effects if the following interactions are allowed: 

qDH d + Le c H d + Dqi + L£e c , (3) 

where H d stands for the MSSM down-type Higgs field. In this case, one has to assume the 
couplings in Eq.Q to be O(10~ 4 ) to obtain the correct neutrino masses. Since our goal in this 
paper is to investigate resonant production of diquark at the LHC, we desire to have at least 
O(O.l) strength couplings of diquarks with the SM quarks. 

Thus, it seems natural to forbid rapid proton decay and large neutrino masses by a sym- 
metry. Clearly, if odd i?-parity is assigned to the vectorlike fields, the couplings in Eq.(j2]) are 
all forbidden, which clearly is not our goal. It was shown in Ref. [9J that a (nB + mL)/2 
£7(1) symmetry, where B and L are the baryon and lepton numbers of the fields and n and m 
are integers, can allow us to retain the "diquark couplings", Dqq and Du c d c , but forbid the 
"leptoquark couplings", Du c e c and Dqi. As a consequence, the rapid nucleon decay operators 
are forbidden. If odd i?-parities are assigned for L + L and even i?-parity for D + D, then all 
the terms in Eq.([3]) can be forbidden and therefore neutrino masses will not be generated from 
the loop effects. 

Let us describe in more detail how the (nB + mL)/2 symmetry works. The charge assign- 
ments for quarks, leptons, Higgs and vector like fields under (nB + mL)/2 symmetry are given in 
Table [H The U(l) symmetry is anomalous and the Green-Schwarz mechanism [16] is applied. 
The SM singlet field S acquires a vacuum expectation value (VEV) close to the string scale M st 
from the Fayet-Illiopoulos D-term potential. If n + m is an odd integer, the AB = AL = ±1 
operators (e.g., qqqi, u c d c u c e c , qqu^e^) are all forbidden, as long as no field with fractional 
charge of (nB + mL)/2 symmetry acquires a VEV. The couplings of quarks and leptons with 
D + D fields are expressed as 

Xjj X jj X jj-\-n-\-m X jj-\-n-\-m 

(■£) T Dm + (£) T bM ° + (£) DuV + (£) Diit (4) 
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If the exponents of S are not integers, the respective couplings are not allowed. As one can see, 
the diquark couplings are allowed only if n + m is an odd number and Xjj is even. 

It is interesting to note that the (nB + mL)/2 symmetry forbids rapid nucleon decay, but 
it can allow baryon number violating processes like the AB = 2 n-n oscillations. In fact, the 
appropriate operators responsible for n-n oscillations are effectively generated [9] if there are 
MSSM singlet fields N : (1, 1) , in addition to the diquark fields. The typical n-n oscillation 
operators are u c d c d c u°d c d c and qqd c ^qqd c \ and the diquark couplings, Du c d c and Dqq, play a 
crucial role in generating them at tree level. Through such couplings, the diquark can decay 
into pairs of up-type and down-type quarks and generate a hadronic dijet or a single top quark 
plus a hadronic jet. 

In order to induce the n-n oscillation operators, one can also add a pair of vectorlike fields 
U : (3, l)-2/3 + U : (3, 1)2/3, an d a diquark coupling \Jd\d C y Since the representations are 
included in 10 + 10 multiplets of SU(5) symmetry, the MSSM gauge coupling unification can 
be maintained if the complete contents of 10 + 10 representations are introduced. Through 
the diquark coupling, the U diquark decay produces only hadronic dijet. Although at high 
energy hadron colliders the dijet signal from the resonant production of the diquarks seems to 
be overwhelmed by QCD background, the large luminosity at the Tevatron allows us to give a 
fairly severe lower bound on the diquark mass [17] . Depending on the couplings, the diquarks 
can be discovered at the LHC through the dijet events [18]. Noting that the D + D diquarks 
can decay into single top quark plus a hadronic jet, we concentrate in this paper on collider 
phenomenology of these diquarks and their decays into top quarks. We find that this process 
involving the top quark gives a lower bound on the diquark mass which is more severe than 
that obtained from the dijet process. Furthermore, the final state top quarks are useful in 
identifying the baryon number of the diquark and a chiral structure of the coupling between 
the diquarks and quarks. 

3 Phenomenological constraint on the diquark couplings 

The magnitude of the coupling between the diquarks and quarks is clearly important for the 
production and decay modes of the diquark fields at the LHC, and we therefore first study the 
phenomenological constraints on these couplings. We define the diquark coupling constants as 

-K^Dq^j + RijDUid'j, (5) 

where i,j are flavor indices. The most important constraint on the coupling constants and the 
vectorlike mass of D, D arise from the oblique corrections for the weak gauge bosons (S and 
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T parameters) [19, 20J. The mass and couplings of the diquarks are also constrained from the 
Tevatron experiment, which will be discussed in the next section. 

In SUSY models, the diquark couplings can induce flavor non-universality in SUSY breaking 
squark mass matrices through the renormalization group (RG) evolutions, even if universality 
is assumed at the unification scale. The flavor violation effects can generate processes via 
the flavor changing neutral current (FCNC) (e.g., meson-antimeson mixings and b — » s(cf)7) 
through loop diagrams involving SUSY particles. Therefore, generic O(0. 1) coupling coefficients 
for the first and second generations are disfavored. In addition to such indirect effects of FCNC, 
the diquark couplings can generate FCNC loop diagrams in which the diquark fields propagate. 

We note that D°-D° meson mixing is generated even at the tree level if a sextet diquark 
field A : (6, 1)| and the coupling Au^u^ are considered. However, in the case of triplet diquarks 
(including U+U), all FCNC processes are generated at the loop level, and the essential features 
for finding the experimental bounds are similar to the MSSM case [21] . Thus the matrix 
elements of kk\ kk, etc., are constrained for a given mass spectrum. 

If the diquark coupling matrix elements are hierarchical, similar to the Yukawa couplings 
for quarks and leptons, the FCNC effects can be adequately suppressed. However, for the 
resonant production of diquarks at the LHC, we prefer a parameter region where the couplings 
for the first generation K\\ (or K12) are of 0(0. 1). We make the following simplifying ansatz for 
minimal flavor violation: 

^ij = fcSij, Rij = R§ij. (6) 

In this case, the RG evolutions do not generate any new sources of FCNC in the squark mass 
matrices. The new direct loop contributions via the diquark fields, however, are generated from 
the left-handed coupling hiDqq because of the Cabibbo-Kobayashi-Maskawa (CKM) mixings. 
The contributions to FCNC processes under the assumption in Eq.(j6]) is similar to the chargino 
loop diagram (except for the charge and couplings), and thus the amount of the contribution is 
not significant compared to that of the SM when flavor universality of the squark squared masses 
is assumed at a unification scale. Therefore, for minimal flavor violation, the experimental 
constraint is satisfied as long as the diquarks are heavier than about 300 GeV, and k < 1. 
Actually, if the diquark coupling is more than O(0. 1), the diquark should be heavier than 
about 300 GeV from the Tevatron result, as we will see in next section. 

Any deviation from Eq.([6]) is severely constrained from K-K mixings, especially for K n - n 2 2, 
K12 and R12 ^ O(10~ 3 ) — O(10~ 4 ). The B^ s -Bd, s mixings and b — > s(d) r y decay give k 23 , R23 
^ O(10~ 2 ), and K13, Rys ^ O(10 -3 ). Such small contributions in the flavor off-diagonal elements 
are not crucial for the collider studies, and so will be disregarded in what follows. 

Through the RG running, the squark masses become smaller at low energy, for a given 
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boundary condition, as the diquark couplings get larger. Under the assumption in Eq.flSJ), the 
SUSY breaking squared scalar masses of diquarks become negative when the coupling is 0(1) 
at low energy due to the fact that the wave function renormalizations of the diquark fields 
are proportional to Tr (kkJ) and Ti^RR^). To avoid color and charge symmetry breaking, an 
appropriate size vectorlike mass of the diquark them needs to be added, in which case the 
fermionic partners of the diquarks are heavier than the scalar diquarks. If the SUSY breaking 
diquark squared masses are negative, the squarks do not become tachyonic even if the diquark 
coupling is large. If the diquark couplings are about 0.3 at low energy, the squared scalar 
mass will not become negative due to the gluino loop contribution and the fermionic diquark 
is lighter than the scalar diquark. 



4 Collider phenomenology 

Since the scalar diquarks (D and D) couple with a pair of quarks, they can be produced in the 
s-channel at the Tevatron and the LHC through annihilation of a pair of quarks. As a signature 
of scalar diquark productions at hadron colliders, we concentrate on its decay channel which 
includes a single top or anti-top quark in the final state. Because of its mass, the top quark 
decays electroweakly before hadronizing. Due to this characteristic feature not shared by the 
other quarks, the top quark can be an ideal tool [22] to probe new physics beyond the standard 
model [23]. 

We consider here the scalar diquark (D) production and the analysis for D is basically the 
same, except for the chirality of the initial and final quarks. The fundamental processes are 
ud —7- D —¥ tb (ud —7- -^-ib for the scalar anti-diquark production). The cross section at the 
parton level is given by 

da(ud -> D -> tb) _ |k| 4 (s - m 2 ) 2 

dcos6 = 48tt (1 - M 2 D ) 2 + M 2 D T 2 tot ' ^ 

Here, we have neglected all quark masses except for top quark mass (m t = 173.1 GeV [21]), 9 
is the scattering angle, and the total decay width r tot of the scalar diquark is the sum of its 
partial decay widths, 

T(D^ud,cs) = ^-\k\ 2 M D , 
on 

r(^tt) = ± W M D (i-g) 2 . (8) 

Note that the cross section is independent of the scattering angle because the diquark is a 
scalar. 
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At the Tevatron, the total production cross section of a quark pair (tb) through the scalar 
diquark in the s-channel is given by 

a(pp -> tb + X) = J dxi J dx 2 J dcos9 [f u (x 1 ,Q 2 )f d -(x 2 ,Q 2 ) + f u (x 1 ,Q 2 )f d (x 2 ,Q' 

da(ud — >■ D — )■ tb; s = xix 2 E^ MS ) 

x d^e ' {} 

where f q denotes the parton distribution function for a quark q in a proton, and Eqms is the 
collider energy. Note that one parton distribution function is for the up (down) quark and the 
other is for the sea down (up) quark. At the Tevatron the production cross sections of D and 

are the same, reflecting the fact that the total baryon number of the initial state pp is zero. 

At the LHC, the total tb production cross section is given by 



o(pp — > tb + X) = J dx\ J dx 2 J dcosO 



x 2 Uxi , QV-fe, Q') d " iUd ^ D ", * 7 XrZ2EU) (10) 

acos# 

Here, both parton distribution functions are for the valence quarks, corresponding to a proton- 
proton system at the LHC. The total production cross section of tb through the scalar anti- 
diquark production in the s-channel is obtained by replacing with the appropriate parton distri- 
bution function for anti-quarks. The initial pp state has a positive baryon number, so that the 
production cross section at the LHC of D is much larger than the one for . The dependence 
of the cross section on the final state invariant mass M inv is given by 

da(pp ->tb + X) _ f 1 4M inv 2 / M s 



dM; 



mv 

B CMS 



f 1 AM- ( M 2 \ 

/ 2 dx— F-Mz, Q 2 )f d — f-, Q 2 a(ud ^D^tb-s = M 2 m ).(ll) 

J ./ nv X - tj CMS \ Xljj CMS J 



Let us first examine the lower bound on the scalar diquark mass from the Tevatron. The 
production of single top quark at the Tevatron at y/s = 1.96 TeV has been observed, and the 
measured cross section is consistent with the standard model prediction [25]. We consider the 
cross section [26] . 

a(pp^tb + X) = 1.05 ±0.81 pb, (12) 

from s-channel W boson, and use it as a constraint on tb production cross section through the 
scalar diquark in the s-channel. Any possible new physics should be in the uncertainty range 
of this cross section, and so we take the bound to be 

a(pp -»■ D ->tb) < 0.81pb. (13) 
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M D /GeV 

Figure 1: The cross section of tb production at the Tevatron with -Ecms = 1-96 TeV mediated 
by the scalar diquark in the s-channel. Solid lines correspond to k — 0.2, 0.3, 0.4 from bottom 
to top. The horizontal dotted line shows the Tevatron bound. 

In our numerical analysis, we employ CTEQ5M [27] for the parton distribution functions with 
the factorization scale Q = m t = 173.1 GeV. Figure 1 shows the total cross section of tb 
production as a function of the scalar diquark mass, with -Ecms = 1-96 TeV. The lower bound 
on the diquark mass is found to be 352, 415, and 459 GeV, for k = 0.2, 0.3, 0.4, respectively. For 
these values of k, we also evaluate the dijet production cross section via the resonant diquark 
production and find that the signal cross section is below the observed upper limit [17] and, 
hence, there no lower bound is obtained. We also find that for a sufficiently large value of k, 
the lower bound on the diquark mass obtained by the observed cross section of single top quark 
production is always more severe than that from the dijet production cross section. 

Next we investigate the scalar diquark and anti-diquark productions at the LHC with 
Eqms — 14 TeV. The differential cross sections for each process for k = 0.3 and with Mp = 600 
GeV and 1 TeV are depicted in Figure 2. Here, we compare the single top quark production 
via the scalar diquark with the it production cross section in the standard model. At the LHC, 
the main background for the single top quark production is it production [28], and it is easy to 
misidentify it events as single top events. So, for a conservative analysis, we compare our signal 
to it production. The single top production cross section is smaller than the it production cross 
section at the LHC by factor of about 7. 

We can see that the peak cross sections for tb production exceeds the standard model cross 
section, while the ib cross section by comparison is lower (however, the ib cross section deviates 
sizably from the standard model cross section). This discrepancy between the production cross 
sections of scalar diquark and anti-diquark at the LHC is direct evidence for non-zero baryon 
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Figure 2: The differential cross sections for tb (solid line), tb (dotted line) versus the invariant 
mass of the final states. The left peak corresponds to Md = 600 GeV and the right one to 
Md = 1 TeV. The dashed line is the standard model ti background. Here k = 0.3 



number of the scalar diquark. Note that one can distinguish the top quark from the anti- 
top quark through their semi-leptonic decays. Counting the number of top quark events and 
anti-top quark events would reveal a non-zero baryon number for the scalar diquark. 

The result for the production of the scalar diquark D is the same when we take k = R and 
replace with D. However, there is a crucial difference in the spin polarization of the final 
state top (anti-top) quark. This is because while D couples with the left-handed top quark, 
D* couples with the right-handed top quark. Since the top quark decays before hadronizing, 
the information of the top quark spin polarization is directly transferred to its decay products 
and results in significant angular correlations between the top quark polarization axis and the 
direction of motion of the decay products [29]. It has been shown that measuring the top 
spin correlations can increase the sensitivity for detecting a new particle at the Tevatron [30J 
and the LHC [31] • For the scalar diquark production, we can distinguish between D and 
productions by measuring the polarization of top quarks from their decays. It is interesting 
that only the right-handed top quark is produced by the scalar diquark decay, while the 
single top quark produced through the standard model process is purely left-handed. 

We have discussed the scalar diquark production in the s-channel, where the production 
cross section is controlled by the diquark coupling. Because they carry color, diquarks can 
be pair produced through QCD processes, and this production is independent of the diquark 
coupling. For the sextet diquark, the pair production at the LHC has been investigated in 
detail [TT]. Although the analysis there focused on tttt final state, a similar analysis here for 
the pair production of triplet scalar diquarks is applicable. The total pair production cross 
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section of the triplet diquark is about an order of magnitude smaller than the sextet diquark 
case because of the difference of the SU(3) group factors (see Fig. 2 in [TT]). 

If kinematically allowed, the scalar diquark can also decay into a pair of squarks through the 
soft SUSY breaking trilinear coupling. In this case, the scalar diquark resonance can enhance 
the production cross section of squarks and can give some advantage to sparticle searches at 
the LHC. 

Fermion diquarks can be pair produced through QCD processes and the production cross 
section is of the same order of magnitude as the production cross section of the scalar diquark 
pair. Since they are i?-parity odd particles, the decay process of fermion diquarks is similar to 
the usual sparticle decay processes. For example, a fermion diquark can decay into squarks and 
quarks, and this process is similar to gluino decay, except for the difference arising from baryon 
number. Since the fermion diquark baryon number is twice that of the quark, it is crucial to 
identify its baryon number from its decay products. However, this can be quite challenging. 

5 Conclusion 

We have investigated a class of SUSY models which contain color triplet diquark fields with 
masses of around a few hundred GeV to one TeV. The couplings between the diquark and a 
pair of quarks play a crucial role in observable n-n oscillation which may be tested in the near 
future. Based on a relatively simple model recently proposed in [9], we have investigated the 
collider phenomenology of the diquark fields. The phenomenological constraints require that 
the couplings between the diquarks and a pair of quarks should be almost flavor-blind and their 
magnitude should be less than unity. For the diquark couplings satisfying these constraints, we 
have analyzed s-channel scalar diquark production and its subsequent decay into the top and 
bottom quarks at hadron colliders. We first investigated the lower bound on the scalar diquark 
mass from the current observation of single top quark production at the Tevatron. We find 
that the diquark mass Mp ^ 415 GeV for a diquark coupling constant of 0.3. At the LHC, the 
differential cross section of the scalar diquark production shows a resonance peak for a diquark 
coupling of 0(0.1). The difference in cross sections for scalar diquark and anti-diquark at the 
LHC, which can be identified by comparing the number of events with top and antitop quarks, 
provides direct evidence of nonzero baryon number of the diquark. The scalar diquarks D and 
are identical except for the chirality of the quarks to which they couple. Measuring the spin 
polarization of top quarks offers a way to distinguish between these diquarks at the LHC. 



11 



Acknowledgments 

This work is supported in part by the DOE Grant No. DE-FG02-91ER40626 (I.G., Y.M. and 
Q.S.) and GNSF Grant No. 07_462_4-270 (I.G.). N.O. would like to thank the Particle Theory 
Group of the University of Delaware for hospitality during his visit. 



References 

[1] See, for instance, E. Witten, Nucl. Phys. B 258, 75 (1985); G. Lazarides, C. Panagio- 
takopoulos and Q. Shafi, Phys. Rev. Lett. 56, 432 (1986); G. Lazarides, P. K. Mohapatra, 
C. Panagiotakopoulos and Q. Shafi, Nucl. Phys. B 323, 614 (1989). 

[2] V. D. Angelopoulos, J. R. Ellis, H. Kowalski, D. V. Nanopoulos, N. D. Tracas and 
F. Zwirner, Nucl. Phys. B 292, 59 (1987); J. L. Hewett and T. G. Rizzo, Phys. Rept. 
183, 193 (1989); S. F. King, S. Moretti and R. Nevzorov, Phys. Rev. D 73, 035009 (2006) 
|arXiv:hep-ph/0510419] ; J. Kang, P. Langacker and B. D. Nelson, Phys. Rev. D 77, 035003 
(2008) jarXiv:0708.270Tl [hep-ph]] . 

[3] K. S. Babu and J. C. Pati, Phys. Lett. B 384, 140 (1996) | arXiv:hep-ph/9606215] ; 
M. Bastero-Gil and B. Brahmachari, Nucl. Phys. B 575, 35 (2000) |arXiv:hep-ph/9907318] ; 
J. L. Chkareuli, I. G. Gogoladze and A. B. Kobakhidze, Phys. Rev. Lett. 80, 912 (1998); 
J. L. Chkareuli, C. D. Froggatt, I. G. Gogoladze and A. B. Kobakhidze, Nucl. Phys. 
B 594, 23 (2001) larXiv:hep-ph/0003007| ; T. Li, D. V. Nanopoulos and J. W. Walker, 
larXiv:0910. 0860 [hep-ph]. 

[4] K. S. Babu, I. Gogoladze, M. U. Rehman and Q. Shafi, Phys. Rev. D 78, 055017 (2008) 
|arXiv:0807.3055l [hep-ph]]: S. P. Martin. laTXTv:0910.2732l [hep-ph]: P. W. Graham, A. Is- 
mail, S. Rajendran and P. Saraswat. larXiv:0910.3 020 [hep-ph]; K. S. Babu, I. Gogoladze 



and C. Kolda, arXiv:hep-ph/0410085 



[5] T. Moroi and Y. Okada, Phys. Lett. B 295, 73 (1992); Mod. Phys. Lett. A 7, 187 (1992). 
[6] S. W. Ham, S. a. Shim and S. K. OH. larXiv: 1001 .11291 [hep-ph], 
[7] T. Ibrahim and P. Nath. larXiv:1001.0231l [hep-ph]. 



[8] Z. Chacko and R. N. Mohapatra, Phys. Rev. D 59, 055004 (1999) |arXiv:hep-ph/9802388|; 
K. S. Babu and R. N. Mohapatra, Phys. Lett. B 518, 269 (2001) [arXiv:hep-ph/0108089] ; 
K. S. Babu, P. S. Bhupal Dev and R. N. Mohapatra, Phys. Rev. D 79, 015017 (2009) 



12 



[9 
[10 



[11 

[12 
[13 

[14 
[15 

[16 

[it; 

[18 

[19 

[20 
[21 



|arXiv:0811.34TTl [hep-ph]]; B. Dutta, Y. Mimura and R. N. Mohapatra, Phys. Rev. Lett. 
96, 061801 (2006) |arXiv:hep-ph/0510291] . 

M. A. Ajaib, I. Gogoladze, Y. Mimura and Q. Shan, Phys. Rev. D 80, 125026 (2009) 
jarXiv:0910.1877l [hep-ph]]. 

E. DelNobile, R. Franceschini, D. Pappadopulo and A. Strumia, Nucl. Phys. B 826, 217 
(2010) [arXiv:0908.1567 [hep-ph]]; T. Han, I. Lewis and T. McElmurry, arXiv:0909.2"666l 
[hep-ph]; A. Arhrib, R. Benbrik and C. H. Chen. larXiv:0911.4875l [hep-ph]. 

C. R. Chen, W. Klemm, V. Rentala and K. Wang, Phys. Rev. D 79, 054002 (2009) 
jarXiv:0811.2T05l [hep-ph]]. 

H. Tanaka and I. Watanabe, Int. J. Mod. Phys. A 7, 2679 (1992). 

R. N. Mohapatra, N. Okada and H. B. Yu, Phys. Rev. D 77, 011701 (2008) 
jarXiv:0709. 14861 [hep-ph]]. 

J. C. Pati and A. Salam, Phys. Rev. D 10, 275 (1974). 

H. Nishino et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 102, 141801 (2009) 
}arXiv:0903. 0676 [hep-ex]]. 



M. B. Green and J. H. Schwarz, Phys. Lett. B 149, 117 (1984); Nucl. Phys. B 255, 93 
(1985); M. B. Green, J. H. Schwarz and P. C. West, Nucl. Phys. B 254, 327 (1985). 

T. Aaltonen et al. [CDF Collaboration], Phys. Rev. D 79, 112002 (2009) [ arXiv:0812.4"036l 
[hep-ex]]. 

S. Atag, O. Cakir and S. Sultansoy, Phys. Rev. D 59, 015008 (1999); O. Cakir and M. Sahin, 
Phys. Rev. D 72, 115011 (2005) |arXiv:hep-ph/0508205] . 

L. Lavoura and J.P Silva, Phys. Rev. D 47, 2046 (1993); N. Maekawa, Phys. Rev. D 52, 
1684 (1995). 

C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 

F. Gabbiani and A. Masiero, Nucl. Phys. B 322, 235 (1989); J. Hagelin, S. Kelley and 
T. Tanaka, Nucl. Phys. B 415, 293 (1994); F. Gabbiani, E. Gabrielli, A. Masiero and 
L. Silvestrini, Nucl. Phys. B 477, 321 (1996) |arXiv:hep-ph/960438 7|. 



[22] C. T. Hill and S. J. Parke, Phys. Rev. D 49, 4454 (1994) [arXiv:hep-ph/9312324| . 



13 



[23] K. Agashe, A. Belyaev, T. Krupovnickas, G. Perez and J. Virzi, Phys. Rev. D 77, 015003 
(2008) |arXiv:hep-ph/0612015] ; B. Lillie, L. Randall and L. T. Wang, JHEP 0709, 074 
(2007) |arXiv:hep-ph/0701166] ; D. Choudhury and D. K. Ghosh, Int. J. Mod. Phys. A 23, 
2579 (2008) |arXiv:0707.2074l [hep-ph]]. 

C. Vellidis [CDF Collaboration]. larXiv:0910.3392i [hep-ex]. 

T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 101, 252001 (2008) 
|arXiv:0809.258Tl [hep-ex]]; V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 103, 
092001 (2009) |arXiv:0903.0850l [hep-ex]]; T. Aaltonen et al. [CDF Collaboration], Phys. 
Rev. Lett. 103, 092002 (2009) [ arXiv:0903.0885l [hep-ex]]. 

C. E. Gerber. laTXiv:0909.4794l [hep-ex]. 

H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12, 375 (2000) 
|arXiv:hep-ph/99032'82] ; 

T. Stelzer, Z. Sullivan and S. Willenbrock, Phys. Rev. D 58, 094021 (1998) 
|arXiv:hep-ph/98073'40| . 

For studies at the LHC see, for example, F. Hubaut, E. Monnier, P. Pralavorio, K. Smolek 
and V. Simak, Eur. Phys. J. C 44S2, 13 (2005) larXiv:hep-ex/0508061]. 



R. M. Harris, C. T. Hill and S. J. Parke, arXiv:hep-ph/9911288, 

M. Arai, N. Okada, K. Smolek and V. Simak, Phys. Rev. D 70, 115015 (2004) 
|arXiv:hep-ph/0409273| ; Phys. Rev. D 75, 095008 (2007) |arXiv:hep-ph/0701155| ; Acta 
Phys. Polon. B 40, 93 (2009) [arXiv: 0804 .37401 [hep-ph]] ; M. Arai, N. Okada and K. Smolek, 
Phys. Rev. D 79, 074019 (2009) [arXiv:0902.04T8l [hep-ph]]. 



14 



